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IV 


1.  Introduction 


In  Mechanics  of  Continua,  Eringen^  presents  a  large  deformation  theory  for  eleetromagneties 
coupled  with  continuum  mechanics.  While  similar  approaches  have  been  taken  by  other  authors 
in  the  literature,  there  is  no  consensus  on  how  to  treat  electromagnetic  fields  in  solids  subjected  to 
finite  deformations.  The  discrepancies  among  these  approaches  almost  always  include 
transformation  rules  for  polarization,  though  there  are  also  inconsistencies  for  the  transformation 
rules  of  the  other  fields  as  well.^^^ 

Eringen  presents  an  application  of  his  theory  in  the  case  of  a  simple  shear  deformation.  The 
boundary  value  problem  (BVP)  solved  corresponds  to  an  infinitely  large  parallel  plate  capacitor 
with  an  isotropic  dielectric  where  one  plate  is  sheared  relative  to  the  other  as  shown  in  the  Eigure. 
Inertial  and  velocity  effects  are  neglected.  This  BVP  is  a  simple,  easily  solved  problem  that  can 
be  used  to  assess  different  hyperelastic  theories,  including  Eringen’s,  since  in  a  hyperelastic 
theory,  the  field  quantities,  e.g.,  polarization,  electric  field,  and  stress,  can  be  determined  in  the 
material  (Lagrangian)  frame,  and  the  spatial  (Eulerian)  frame  as  a  function  of  deformation  once  a 
form  of  the  free  energy  function  has  been  assumed. 

In  this  report,  we  detail  the  solution  of  the  simple  shear  example  problem  using  the  hyperelastic 
theories  of  Eringen'  and  Clayton.^  Both  theories  predict  a  polarization  response  that  is  not 
collinear  with  the  electric  field  in  the  spatial  frame.  Eurthermore,  this  effect  is  first  order  in  the 
deformation  and  does  not  reduce  to  accepted  small  deformation  theory  where  there  is  no 
distinction  between  the  spatial  and  material  frames. 


y,Y 


Figure.  Schematic  of  the  dielectric  slab 
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2.  Nonlinear  Theory  and  Material  Frame  Representation 

In  general,  the  deformation  at  a  material  point  (see,  for  example,  Malvern^)  is  defined  by  a 
transformation  from  the  material  coordinates,  X^,  to  the  spatial  coordinates,  x\  through  the 
mapping 

=  x\X^)  .  (1) 

Throughout  this  report,  lowercase  subscripts  refer  to  the  spatial  system  while  uppercase 
subscripts  refer  to  the  material  system.  The  mapping  x\X^)  is  assumed  to  be  bijective  so  that  the 
inverse  relationship  X^  =  X\x^)  exists.  The  mapping  is  also  assumed  to  be  smooth  and 
differentiable.  The  deformation  gradient  tensor,  x^j  =  is  defined  through 

dx^  =  x^'jdX^  ,  (2) 

with  the  inverse  deformation  gradient  tensor  given  by  X{.  The  Green’s  strain  tensor  is  defined  in 
terms  of  the  deformation  gradient  tensor  by 

Cij  =  XjQijX^j  ,  (3) 

where  gij  is  the  metric  tensor  of  the  spatial  frame. 

2.1  Electrostatics  and  Boundary  Conditions 

For  this  example,  electrostatics  is  assumed,  and  all  the  governing  equations  are  presented  with 
respect  to  the  spatial  reference  frame.  The  pertinent  Maxwell’s  equations  are  in  either  vector  or 
indicial  notation  where  superscripts  and  subscripts  refer  to  contravariant  and  covariant 
components,  respectively,  given  by 

VxE  =  0,  (5) 

V\S\ 

where  the  standard  rule  of  summation  over  repeated  indices  applies.  The  relationship  between  the 
electric  displacement,  Di,  electric  field,  Ei,  and  polarization.  Pi,  is 

Di  =  e^Ei  +  Pi  ,  (6) 


2 


where  eo  is  the  permittivity  of  free  space.  In  the  Heaviside-Lorentz  units  used  by  Eringen,  eo  is 
unity. 

The  boundary  conditions  are  determined  from  the  jump  conditions  on  the  material  interfaces 

n-[B]=Wf  ,  (7) 

n  X  [E]  =  0  ,  (8) 

where  Wf  is  the  surface  charge  and  n  is  the  normal  to  the  surface.  Following  Eringen,  a  constant 
uniform  electric  field  is  applied,  which  satisfies  Eq.  5.  This  reduces  Eq.  4  to 

V  ■  P  =  0  .  (9) 

Consequently,  the  polarization  is  spatially  constant  as  well.  If  the  applied  electric  field,  in  the 
spatial  frame  as  the  material  is  deforming,  only  has  an  component,  the  components  of  the 
electric  field  in  the  spatial  frame  are 

E,  =  0,  Ey  =  E,  =  0  ,  (10) 

where  x  =  x^,y  =  x'^,  and  z  =  x^  have  been  used  for  notational  convenience. 

2.2  Mechanical  Deformation  and  Boundary  Conditions 

As  mentioned,  the  region  of  interest  is  defined  to  be  an  infinite  isotropic  homogeneous  dielectric 
slab  subjected  to  a  simple  shear  deformation  defined  by 

+  kX^  , 

x‘^  =  X^,  (11) 

x^  =  X\ 

where  k  is  the  amount  of  shearing.  The  associated  deformation  gradient  tensor,  x\,  and  Jacobian, 
J,  are 

'l  k  o' 

[</]=  0  10,  (12) 

_0  0  1_ 

J  =  det{xij)  =  1  .  (13) 


3 


3.  Eringen’s  Hyperelastic  Theory 


Eringen^  presents  a  general  nonlinear  theory  for  isotropic,  elastic  dielectrics,  including  both 
material  and  geometric  nonlinearities 7  In  this  section,  this  hyperelastic  theory  is  presented  and 
applied  to  the  scenario  of  simple  shear  of  a  dielectric  material,  which  was  described  in  the 
previous  section.  A  form  of  the  free  energy  is  assumed  such  that  the  material  is  linear  and 
isotropic  in  the  material  frame.  This  assumed  form,  coupled  with  simple  shear  deformation,  is 
then  used  to  determine  the  finite  deformation  electric  field  in  the  spatial  frame.  Eringen’s  theory 
was  developed  assuming  both  the  spatial  and  material  frames  are  Cartesian.  Therefore,  all  the 
superscripted  indices  in  the  preceding  sections  can  be  lowered  to  subscripts.  Also,  in  Cartesian 
frames,  the  metric  tensors  are  the  identity  tensor,  e.g.,  gij  =  5ij. 

The  deformation  gradient  tensor  is  used  by  Eringen^  to  transform  the  mechanical  and 
electromagnetic  fields  between  the  spatial  and  material  reference  frames.  Eor  electrostatics, 
Eringen  defines  the  material  to  spatial  frame  transformations  as 


Ek  —  ExXx^k  , 

Pk  =  -fP-KXk,K  , 


(14) 


where  Ek  and  Pk  are  the  electric  field  and  polarization  represented  in  the  material  frame.  J  is  the 
determinant  of  the  deformation  gradient  tensor. 


3.1  Thermodynamics 


According  to  general  thermodynamic  theory,  the  material  properties,  stress  state,  and  polarization 
can  be  determined  from  the  deformation  and  the  free  energy,  4/,  which  is  the  energy  available  in 
the  system  for  mechanical  work.  Eor  electrostatics,  the  free  energy  is  assumed  to  only  be  a 
function  of  the  Green’s  strain,  Cjj,  and  the  material  electric  field,  Ej. 


Po^  =  J:{Cij,Ej)  .  (15) 

Eor  an  isotropic  material,  S  must  be  objective  to  rotations  and  therefore  can  be  written  in  terms  of 
the  principal  traces.  (Eor  a  more  in-depth  treatment  of  objectivity  requirements  for  isotropic 
functions,  see  Itskov.^) 
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The  polarization  can  be  obtained  in  terms  of  the  free  energy  by 


-  2^CuEj  -  2^CijCjkEk  , 

(7/4  <^-^6  ^^8 


(16) 


where 

J4  =  E  ■  E,  Je  =  E  ■  C  ■  E,  Jg  =  E  ■  ■  E  .  (17) 

Since  it  is  not  the  focus  of  the  current  work,  we  simply  mention  that  the  stress  can  also  be 
determined  in  general  form  as  a  function  of  the  energy  if  needed.^ 


3.2  Isotropic  Linear  Dielectric  Material 


Assuming  the  material  is  an  isotropic,  linear  elastic,  homogeneous  dielectric,  the  form  of  S, 
simplified  from  Eringen,^  is 

E  =  2^^^^  2^(,l2  +  -OiglA  , 


(18) 


where  the  are  material  properties,  and  Ji,  I2,  and  I4  are  the  principal  traces 

Ii  =  tr{C),  l2  =  tr{C^),  h  =  E-E.  (19) 

Substituting  Eq.  18  into  Eq.  16  gives 


11/  —  —a^Ej  .  (20) 

The  material  parameter  —as  is  related  to  the  dielectric  constant.  The  parameters  ai  and  can  be 
shown  to  be  Eame’s  constants.  With  the  form  of  the  energy  given  by  Eq.  18,  the  polarization  does 
not  depend  on  the  strain,  and  the  stress  does  not  depend  on  the  electric  field.  Therefore,  the 
mechanical  and  the  electrostatic  responses  are  uncoupled.  The  polarization  is  also  collinear  with 
the  electric  field  in  the  material  frame. 

Based  on  the  transformations  in  Eq.  14  and  the  boundary  conditions,  the  electric  field  in  the 
material  frame  is 

Ex  =  0,  Ey  =  Ey,  Ez  =  0  .  (21) 

Since  the  electric  field  is  known  from  Eq.  21,  the  polarization  in  the  material  frame  is 

Ex  =  0,  Ily  =  —agEy  =  —agEy,  Ilz  =  0  .  (22) 

Erom  Eq.  14,  the  polarization  components  expressed  with  respect  to  the  spatial  frame  can  be 
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determined  from 


’l 

k 

o' 

f  0] 

P].= 

0 

1 

0 

< 

nv 

_0 

0 

1_ 

io  J 

The  spatial  components  of  the  polarization  are  thus 


P. 

Py 

Pz 


(  kUy 

I  ny 

I  0 


1  -"f  1 


(23) 


(24) 


It  is  seen  that  the  final  form  for  the  polarization  in  the  spatial  frame  is  not  collinear  with  the 
electric  field  due  to  the  Px  component  depending  on  the  magnitude  of  the  shear  deformation 
through  k. 


3.3  Discussion  of  Eringen’s  Theory 


The  finite  deformation  electrostatic  theory  presented  by  Eringen  leads  to  anisotropy  of  the 
dielectric  material  even  when  an  isotropic,  linear,  homogeneous  material  is  assumed.  This 
anisotropy  is  why  the  polarization  and  electric  field  are  not  collinear.  Rewriting  Eq.  24  so  that  the 
effective  material  properties  appear  as  a  matrix. 


-a8(l  +  k‘^) 

—ask 

0 

>  = 

—ask 

—as 

0 

0 

0 

—as 

E?. 


and  defining 


-08(1  +  k"^) 

—agk 

0 


—ask  0 
-as  0 
0  —as 


(25) 


(26) 


it  can  be  seen  that  the  [aju  term  has  a  nonlinear  dependence  on  the  deformation  that  has  no  effect 
when  the  electric  field  is  1 -dimensional  in  the  x'^  direction. 


Also,  note  that  while  the  magnitude  of  the  polarization  in  the  direction  increases  with  the 
deformation,  the  energy  associated  with  the  electric  field  is  constant.  This  is  because  the  energy 
associated  with  the  polarization  and  the  electric  field  is  E  ■  P,  and  since  there  is  no  Ex  component 
of  the  electric  field,  the  Px  component  of  the  polarization  does  not  contribute  to  the  total  energy. 
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4.  Clayton’s  Hyperelastic  Theory 


A  hyperelastic  approach  is  also  presented  for  electroactive  materials  by  Clayton.^’^  Again, 
electrostatics  is  assumed  as  in  Section  2.1,  and  the  applied  deformation  as  defined  in  Section  2.2, 
respectively.  Clayton  developed  the  theory  using  general  curvilinear  coordinates,  so  this  section 
uses  raised  indices  for  contravariant  terms  and  lowered  indices  for  the  covariant  terms. 

The  mappings  between  the  spatial  and  material  fields  defined  by  Clayton  are 


Ea  =  XyEa  , 

Pa  =  xyPa , 


(27) 


where  hats  are  used  for  consistency  with  Clayton.  The  first  of  Eq.  27  is  consistent  with  the 
transformation  Eringen  used  for  the  electric  field,  while  the  polarization  transform  is  different. 
Also,  the  transformation  for  the  electric  field  is  consistent  with  the  electrostatic  version  of 
Earaday’s  law  (Eq.  5)  being  form  invariant  between  the  spatial  and  material  frames.  These 
transformations  are  for  the  covariant  components  of  the  fields,  noted  by  the  position  of  the  indices. 


4.1  Thermodynamics 


The  constitutive  assumptions  are 


i;  =  i;{EAB,pA,0,0,A,^,GA)  , 

E^  =  E^{Eab,  Pa,  0, 0,A,  X,  Ga)  ,  (28) 

TAB  ^  ^ 

where  ip  is  the  free  energy,  E"^  is  the  electric  field,  and  is  the  second  Piola-Kirchhoff  stress 
tensor  (not  symmetric).  Eab  is  the  Eagrangian  strain  defined  as 

Eab  =  -^{Cab  —  Sab)  ■  (29) 

These  assumptions  differ  from  those  of  Eringen,  Eq.  15,  in  that  polarization  is  taken  as  an 
independent  variable,  and  the  electric  field  is  a  dependent  variable,  while  Eringen  assumed  the 
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opposite.  The  relationships  between  the  free  energy  and  the  other  dependent  variables  are 


Ea  —  CabP 
=  po 


dip 


OPb 

dip 


dEAB 


+  JC-^^^PcC-^^^Ed  ■ 


(30) 


As  mentioned  earlier,  we  only  indicate  that  the  stress,  can  be  determined  from  the  free 
energy,  but  we  do  not  consider  this  further.  To  make  the  notation  similar  to  Eringen,  let 


poip  =  T.  . 


(31) 


Now  assume  that 

S  =  +  2®8-^4  ;  (32) 

where 

h  =  fr(E),  h  =  tr(E2),  ^  =  P  ■  P  .  (33) 

This  is  equivalent  to  the  form  of  the  free  energy  for  the  linear  isotropic  material  used  by  Eringen, 
Eq.  18,  except  the  terms  containing  the  electric  field  are  switched  to  polarization,  which  are 
denoted  by  the  use  of  a  tilde  over  the  variable  name.  Note  that  ds  is  related  to  the  dielectric 
properties  but  is  not  equal  to  ag. 

Using  Eq.  32,  the  electric  field  is 

Ea  =  EaAE^  =  EaAE^-^^  ,  (34) 

Po  dPB 

where 

E^  =  x%  EaA  =  E\g,a  .  (35) 

The  EaA  term  appears  in  Eq.  34  because  of  the  need  to  lower  the  contravariant  and  to  perform 
the  pull-back  to  the  material  frame.  Using  the  assumed  form  of  the  free  energy,  Eq.  32,  the 
electric  field,  Eq.  34,  becomes 

1  ..  dy 

Ea  =  -Cab{‘2P^^)  ,  (36) 

with  the  right  Cauchy-Green  strain  defined  in  curvilinear  coordinates  as 

Cab  =  xygabx%  .  (37) 
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From  Eq.  32, 


(38) 


<9S 

— ^  —  ~^8  1 

dh  ^ 

SO  the  electric  field  can  finally  be  determined  by  substituting  Eq.  38  into  Eq.  36, 

Ea  =  -jO^sCaeP^  ■  (39) 

This  is  the  constitutive  relation  between  Ea  and  in  the  material  frame.  Eirst,  the  position  of 
the  indices  on  the  fields  should  be  noted  since  the  covariant  and  contravariant  nature  of  the  fields 
is  important  with  regards  to  which  bases  the  fields  are  referred.  The  implications  for  this  become 
apparent  when  mappings  between  the  frames  are  considered  in  the  next  section.  Second,  Eq.  39 
indicates  that  the  deformation  is  explicitly  included  in  the  constitutive  relationship  between  the 
electric  field  and  polarization  in  the  material  frame. 

4.2  Simple  Shear  Example 

As  before,  a  simple  shear  deformation  under  a  constant,  uniform  electric  field  is  examined  to  see 
how  the  theory  is  applied.  The  deformation  is  as  defined  in  Section  2.2.  Additional  assumptions 
include  that  both  the  spatial  and  material  frames  are  Cartesian,  as  was  also  the  case  with  Eringen’s 
theory,  and  consequently  the  metrics  are 

gab  =  5ab,  =  Gab  =  Sab,  =  5^^  .  (40) 


This  has  the  effect  of  making  the  right  Cauchy-Green  strain.  Gab,  the  same  as  Eringen’s.  Even 
though  Cartesian  systems  are  assumed,  curvilinear  notation  will  still  be  employed  since  it  is 
essential  in  the  discussion  of  the  relationships  between  the  fields. 


We  again  assume  that  the  spatial  electric  field  is  in  the  -direction,  so 


Using  the  transformation  for  the  electric  field,  Eq.  27,  the  material  field  is 


’l 

0 

o' 

[A]  = 

k 

1 

0 

< 

4 

_0 

0 

1_ 

.  0  J 

(41) 


(42) 
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From  Eq.  39, 


(43) 


1 

as 


B 


where  C  is  the  inverse  of  Cab,  which  for  simple  shear  is 


-lABi 


l  +  k^ 

-k  1 
0  0 


k  0 
0 


(44) 


This  gives  a  polarization  in  the  material  frame  of 


—  kEy 
0 


(45) 


The  transformations  between  the  frames  given  by  Eq.  27  are  covariant  in  nature,^  requiring  the 
polarization  to  be  shifted  from  contravariant  components  to  covariant  components  using 


Pa  =  GabP^  =  SabP^  , 


yB 


(46) 


SO  that  the  covariant  components  are  equal  to  the  contravariant  ones, 


-kEy 
0 


(47) 


Using  Eq.  27,  the  spatial  polarization  field  is  then  given  by 


[^■J  = 


-kE° 

^<j(i  +  F)£; 

0 


(48) 


4.3  Discussion  of  Clayton’s  Theory 

Similarly  to  Eringen’s  theory,  the  polarization  is  not  collinear  with  the  applied  electric  field  in 
either  the  material  or  spatial  frames.  In  contrast  with  Eringen’s  theory,  the  -component  of  the 
polarization  has  a  second-order  term  associated  with  the  deformation,  Py  =  apl  +  kpEy. 
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Rewriting  Eq.  48  similarly  to  Eq.  26  so  that  the  dielectric  material  properties  appear  as  a  matrix, 


and  defining 


—  ^k 

OLS 

0 


0 

0 

J_ 

as- 


Ey 


as 


0 


—  ^k 

as 


0 

0 

_j_ 

08- 


the  full  dependency  of  the  polarization  on  the  deformation  can  be  seen. 


(49) 


(50) 


5.  Conclusions 


The  problem  solved  is  analogous  to  an  infinite  dielectric  slab  polarized  by  a  surface  charge  (or 
applied  electric  field).  The  boundary  conditions  for  the  electrostatics  are,  consequently, 
independent  of  the  deformation.  The  distribution  of  surface  charges  relative  to  one  another  is  not 
affected  by  the  shear  deformation.  In  other  words,  a  positive  charge  on  the  upper  surface  of  the 
slab  has  no  preference  for  interacting  with  a  specific  negative  charge  on  the  bottom  surface  of  the 
plate.  There  are  3  issues  that  will  be  discussed  in  relation  to  the  solution  to  the  problem:  the 
anisotropy,  the  lack  of  consistency  between  the  two  theories,  and  the  lack  of  consistency  with 
small  deformation  theory. 

It  is  known  that  for  the  elastic  tangent  stiffness,  under  isotropic,  hyperelastic  assumptions,  a 
material  can  exhibit  deformation-induced  anisotropy. There  are  thermodynamic  constraints  on 
the  free  energy  that  determine  when  the  stiffness  tensor  retains  isotropy  under  deformation.  These 
constraints  limit  how  the  free  energy  can  depend  on  the  deformation.  Similar  thermodynamic 
constraints  are  unknown  for  the  dielectric  material  properties,  but  the  linear  assumptions  made  in 
Eqs.  18  and  32  are  consistent  with  the  constraints  for  isotropy  of  the  stiffness  tensor.  Eurther  study 
must  be  done  to  determine  if  the  dielectric  properties  should  remain  isotropic  under  deformation. 

Second,  the  2  hyperelastic  theories  for  electroactive  materials  give  different  answers  for  the  same 
simple  shear  problem.  While  both  theories  predict  the  development  of  deformation-induced 
anisotropy  (the  [d]i2  and  [q;]2i  terms  in  Eqs.  26  and  50),  nonlinear  effects  appear  in  different 
places.  Eor  Eringen’s  theory,  the  an  term  has  the  nonlinear  dependence  on  deformation  (1  +  k'^). 
The  same  dependency  appears  in  Clayton’s  theory  but  in  the  0:22  term. 
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Finally,  under  the  assumption  of  small  deformations,  the  material  and  spatial  frames  are  the  same. 
For  simple  shear,  this  condition  applies  when  A;  1  in  Eq.  11,  so  that  — )■  0.  Both  theories 
converge  to  the  same  solution  under  this  condition,  with  the  nonlinear  terms  disappearing.  The 
deformation-induced  anisotropy  is  retained,  since  it  depends  linearly  on  k.  This  is  problematic, 
though,  since  there  is  no  difference  between  the  material  and  spatial  frame.  Under  small 
deformations,  an  isotropic  dielectric  material  should  retain  isotropy  under  deformation. 

The  concerns  with  the  solutions  discussed  previously  demonstrate  that  the  theory  for  electroactive 
materials  coupled  with  solid  mechanics  is  not  settled.  One  source  of  the  ambiguity  is  the 
transformations  used  for  mapping  the  electric  field  and  polarization  between  frames.  The 
motivation  for  the  transformation  for  the  electric  field  is  that  Maxwell’s  equations  have  the  same 
form  in  both  the  spatial  and  material  reference  frames.  Other  authors  make  similar  assumptions. 
Since  the  polarization  does  not  explicitly  appear  in  Maxwell’s  equations  (appearing  only  through 
its  relationship  to  E  and  D),  there  is  ambiguity  in  the  literature  over  how  the  polarization  maps 
between  the  spatial  and  material  frames.  The  problem  is  more  fundamental  than  determining  how 
the  polarization  maps  between  frames,  though.  Maxwell’s  equations  are  not  form  invariant  under 
the  types  of  transformations  defined  between  the  spatial  and  material  frames.  This  was  shown  by 
Einstein^^  in  his  theory  of  special  relativity.  To  find  a  fully  consistent  set  of  electromagnetic  field 
transformations,  special  relativity  must  be  invoked  and  then  simplifying  assumptions  made  for 
low-velocity  cases  (e.g.,  see  Wiele  et  al.^^).  More  work  must  be  done  to  develop  a  fully  coupled 
electromechanical  theory  that  addresses  the  inconsistencies  existing  in  the  current  literature. 
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